Design of magnetic textures of ad- nano corrals with an extra adatom 
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It is shown that in antiferromagnetic open or closed corrals of magnetic adatoms grown on surfaces, 
the attachment of a single extra adatom anywhere in the corral impacts on the geometrical topology 
of the nanosystem and generates complex magnetic structures when a magnetic field is applied or a 
magnetic coupling to a ferromagnetic substrate exists. The spin configuration of the corral can be 
tuned to a non-planar state or a planar non-collinear or ferrimagnetic state by adjusting its number 
of sites, the location of the extra adatom or the strength of the coupling to the ferromagnetic 
substrate. This shows the possibility to generate non-trivial magnetic textures with atom by atom 
engineering anywhere in the corral and not only at the edges. 

PACS numbers: 75.75.-c Magnetic Properties of Nanostructures, 75.50.Ee Antiferromagnetics, 75.30.Hx 
Magnetic impurity Interactions. 
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In recent years artificial engineering of molecular 
nanomagnets, magnetic clusters and arrays of magnetic 
adatoms adsorbed on surfaces has emerged for the con- 
struction of entities with rich magnetic properties that 
can be constituents of nanospintronics devices 
These entities can be fabricated directly on surfaces in a 
bottom-up fashion with Scanning Tunneling Microscopy 
(STM) |9f0, and STM is also used to directly measure 
their magnetic properties [l2|, [l3|. The rich magnetic 
properties originate in the exchange couplings between 
the individual magnetic moments [M-Tt] and are of great 
interest for concepts like spin-transfer torque [18-20.] and 
spin chirality [21| on the nanoscale, as well as for po- 
tential applications in quantum computing [25-^. The 
STM measurements allow precise access to the properties 
of the individual magnetic moments, which when com- 
bined with tailored construction of the coupling between 
the spins provide systems where nanomagnets with de- 
sired properties can be synthesized [25|, l26| . In particu- 
lar, arrays of a small number of magnetic atoms on non- 
magnetic metallic substrates are very promising candi- 
dates in this direction. This bottom-up approach is one 
of the main focus research areas of nanoscience where 
various atomic structures, e.g. corrals of adatoms and 
nanowires, are engineered atom by atom. 

It has already been shown theoretically using Density 
Functional Theory (DFT) that adatom nanochains can 
support non-collinear magnetic structures when the ex- 
change interaction between the adatoms is antiferromag- 
netic (AFM) and a coupling to a ferromagnetic substrate 
exists [l^]. In the peculiar case of Mn chains on a Ni(lOO) 
substrate, the AFM coupling within the chain competes 
with the ferromagnetic coupling of the chain to the sub- 
strate leading to an even-odd effect, where the magnetic 
texture depends crucially on the parity of the number of 
atoms in the chain. The Mn chain is in an AFM config- 
uration but if the weak ferromagnetic interaction to the 
substrate is switched on, the spins in the odd numbered 



chain retain their collinearity while the uncompensated 
moment of the chain aligns with that of the substrate. 
Even numbered chains, however, develop a more com- 
plex non-collinear ground state. This even-odd effect has 
been observed recently by STM in short Mn wires on a 
Ni(llO) substrate [llj. 

For such an effect a ferromagnetic substrate is required 
when the exchange interaction is nearest neighbor, there- 
fore strong. However, when the exchange interaction 
is non- nearest neighbor, therefore weaker, for example 
when it is mediated by Ruderman-Kittel-Kasuya-Yosida 
(RKKY) interactions [29'-'3l], an external magnetic field 
can provide the necessary energy to compete with the 
exchange energy. Measurements with STM can be per- 
formed for both kinds of situations: Magnetic nanostruc- 
tures on ferromagnets or on non-magnetic substrates. It 
is pointed out that if the magnetic anisotropy energy is 
strong (a few meV), the RKKY interactions would not be 
able to create non-collinear structures although there is 
competition between the involved magnetic interactions. 

In this paper, AFM wires or corrals with a small num- 
ber of magnetic adatoms are considered, where an extra 
adatom is attached anywhere along the nanostructure. 
It is shown that one more or less adatom is crucial for 
the magnetic properties of such nanostructures, and can 
generate planar or even non-coplanar magnetic configu- 
rations. This opens up the possibility of generating new 
complex magnetic structures whose accessibility depends 
on the precise location of the extra adatom, which is 
not necessarily attached at the edges of the corral, the 
total number of sites of the corral and the strength of 
the coupling to a ferromagnetic substrate or an exter- 
nal magnetic field. Besides the parity of the number of 
adatoms in the corrals, the parity of the site to which the 
extra atom is attached determines not only the magnetic 
texture of the whole nanostructure (open or closed), but 
also the magnetic behavior of the edge atoms for the open 
corrals. This is of crucial importance for building logic 
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gates made up of a few adatoms (see e.g. Q). In the 
following, we will refer to wires which are not necessarily 
straight as open corrals. 

The magnetic properties of the nanostructures are 
modeled with the AFM Heisenberg model (AHM). It is 
mainly considered at the classical level, but it is also 
shown that the classical results compare very well with 
results for typical quantum values of the magnetic mo- 
ments. Previous calculations based on DFT have shown 
that the AFM Heisenberg model is reliable to predict the 
complex magnetic texture in such nanostructures p7j . 
The AHM is considered for corrals of adatoms that form 
open or closed corrals, with an extra adatom attached 
to any location other than the edges. The coupling be- 
tween the individual spins 5^ is J > 0. The interaction 
with a ferromagnetic substrate or with an external mag- 
netic field (in the case of a non-magnetic substrate) has 
strength /i, and the Hamiltonian is: 

AT-l N 
H = J{^Si- Si^i + SL • Se) - s\ + 5|) (1) 

i=l i=l 

where the number of atoms of the corral is A/", with 
L = 2, A/" — 1 the corral adatom away from the edges 
that couples to the extra adatom se- For closed cor- 
rals the term Js\ • must be added to ([T]) and L can 
be any site as there are no open ends and the system is 
translationally invariant without the extra adatom. The 
magnetic field h lies along the z axis and tends to make 
all spins parallel, competing with the AFM configuration 
supported by J. To simplify the discussion the magnetic 
interaction between the extra adatom E and the adatom 
at L is taken to be J, but this is not necessarily the 
case. The result of varying this interaction is discussed 
briefiy in the Supplemental Material. For the case of 
RKKY interactions, small magnetic anisotropy energies 
are required in order to observe a rich magnetic phase di- 
agram, otherwise the adatom moments would be pinned 
to a collinear behavior. This is the case for Fe adatoms on 
a Cu(lll) surface [26]. It is also demonstrated that small 
magnetic anisotropy energies, not considered in the afore- 
mentioned Hamiltonian, lead to negligible effects on the 
magnetic behavior of the nano-corrals (see Supp. Mat.). 
Thus adatoms with low magnetic anisotropy energy, e.g. 
Cr or Mn, are proposed, in contrast to Co or Fe, to be 
deposited on different non-magnetic substrates. Cr and 
Mn have weak magnetic anisotropy so that the AFM 
Heisenberg model without an extra anisotropy term ac- 
curately describes the magnetism of the system. 

DFT calculations predict that for Cr, Mn and Fe de- 
posited on Cu(lll) the magnetic moments are respec- 
tively 4.1, 4.3 and 3.2 \ib [32], while Cr and Mn adatoms 
on a Ni(OOl) or a Fe(OOl) substrate have a magnetic mo- 
ment of the order of 3.5 to 4 /i^, which decreases when 
forming chains due to hybridization [14]. Thus, typically 
Si = 3/2 or 2. It is sufficient to consider the Si as clas- 
sical unit vectors, and the lowest energy configuration is 
found for any h [H, [sj- Typical quantum values are 



also considered [35|, [s^ . It has already been shown that 
for open chains with no extra spins attached the classi- 
cal predictions survive for relatively low values of Si [33] • 
Here similar conclusions are drawn for the dependence of 
the magnetic properties on the magnitude of Si (see also 
Suppl. Mat.). 

For open corrals the exact location of the extra adatom 
is important, due to the lack of translational symmetry 
of the open corral. The magnetic behavior depends on 
the parity of N ^ and also on the parity of the linking 
location L of the extra adatom se- For closed corrals the 
exact linking point is unimportant due to translational 
symmetry, and the magnetic behavior depends only on 
the parity of N . 

Open even corrals. An open even corral with no extra 
adatom attached has nearest neighbor spins pointing in 
opposite directions in the absence of an external field. 
Irrespectively of the point of attachment of the extra 
adatom, the extra spin is not compensated and the total 
magnetization at zero field Mh=o = ^i, corresponding to 
a ferrimagnetic (FI) configuration (Fig. [II[a)). It takes a 
finite magnetic field he to change the FI configuration, up 
to which the only gain in energy comes from the coupling 
to the field. The susceptibility x is discontinuous as he 
is crossed and for higher fields the spins are in a planar 
non-collinear (NC) configuration. The open even corral 
with an adatom attached to it is thus similar to an open 
odd chain [37]. he decreases on the average with N (Fig. 
[21(a)), as Mh=o decreases with respect to the saturation 
magnetization Msat = (A" + l)si with N, = 
he also depends on the parity of L. As L changes from 
even to odd, the length of the largest sub-chain of the cor- 
ral changes its parity from even (which when isolated has 
no FI lowest field configuration) to odd (which possesses a 
FI low field configuration), and correspondingly its zero- 
field magnetization from to Si. It is thus reminiscent 
of an isolated even or odd chain respectively, especially 
since the corral feels the infiuence of the extra adatom 
more strongly around its linking site L, and its infiuence 
gradually weakens when going away from it. Therefore 
the length of the largest odd sub-chain fiuctuates signif- 
icantly with L when the linking point is away from the 
middle, generating the pronounced non-monotonic de- 
pendence of he on the parity of L (Fig. Efa)), due to the 
monotonic dependence of the critical field on the length 
for an isolated (odd) chain. When the extra adatom ap- 
proaches the center of the chain the length of the largest 
odd sub-chain does not change significantly and he tends 
to a constant. 

Open odd corrals. An open corral with an odd num- 
ber of atoms has a single uncompensated spin. When 
an extra adatom is attached Mh=o depends on the par- 
ity of the linking point L. If L is even the extra spin 
is parallel to the uncompensated spin of the open corral 
and Mh=o = 2^^ (Fig. [ijbl)). In contrast, if L is odd 
the extra spin balances out the uncompensated spin of 
the corral and Mh=o = (Fig. [Hb2)). For even L the 
lowest energy configuration is FI, similarly to an open 
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FIG. 1: Lowest energy spin configurations in the absence of a 
magnetic field (red arrows denote spins in the corral, and the 
blue arrow the attached spin), (a) Open even corral with a 
number of adatoms N —Vd and an extra adatom at location 
L — 2. The zero field magnetization is Mh^o = Si. (bl) Open 
odd corral with N = 9 and an even linking point L = 2, 
Mh=o = 2si. (b2) Open odd corral with = 9 and an odd 
linking point L = 3, Mh^o = 0. (c) Closed even corral with 
N = 10, Mh=o = Si. (d) Closed odd corral with AT = 9, 
Mh^o = Si. As in (c) the value of L does not make any 
difference due to translational invariance. 



even corral with an extra adatom. However there are 
now two uncompensated spins, consequently he is big- 
ger (Fig. Efb)). For the same reason he is also bigger 
compared to an open even corral, which also has only a 
single uncompensated spin. For higher fields the lowest 
energy configuration is also NC as in the open even corral 
case (Fig. [3]). The change between the FI and the NC 
configuration also generates a discontinuity in the mag- 
netic susceptibility x (upper left inset of Fig. [3]), with 
nearest-neighbor spins pointing in opposite directions in 
the azimuthal plane for higher fields (lower right inset of 
Fig. [3]). On the other hand, when L is odd, the lowest 
energy configuration is AFM for = and M immedi- 
ately responds to an external field with no susceptibility 
discontinuity (Fig. El^b)), leading to a NC configuration, 
similarly to an even corral with no extra adatoms. Con- 
sequently there are changes in the lowest energy config- 
uration as a function of h only for even L, where the 
size of the largest odd sub-chain of the isolated corral 
decreases with L. Thus the critical field for increasing 
even L increases as the zero-field sub-chain magnetiza- 
tion increases with respect to the sub-chains saturation 
magnetization, in contrast to the non-monotonic result of 
the open even corrals. This shows the importance of the 
exact location of attachment of the extra adatom for the 
magnetic configuration of the corral, leading to an even- 
odd linking effect for open odd corrals, which is absent 
for open even corrals. These clear-cut differences between 
open-even and open-odd corrals provide the opportunity 
for a variety of magnetic structures as the length of the 
corral is varied and the extra adatom is moved along the 
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FIG. 2: Critical fields he scaled to the saturation field hsat for 
which transitions between different configurations take place 
for the corrals, (a): open corrals with even number of adatoms 
N, (b): open corrals with odd N. The horizontal axis shows 
the location L where the extra adatom is attached. The values 



of ■ 



are symmetric with respect to the center of the corral. 



so for an even corral the values for L = 2, . . . , ^ are shown, 
while for an odd corral the values for L = 2, . . . , ^^^^ . A single 
corral is represented by symbols of the same color, and the 
direction where N increases is shown with an arrow. For fixed 
L, decreases with increasing N. For odd corrals there 

' hsat ^ 

is a transition only for even values of L. A clear difference 
in the behavior of he between even and odd corrals is seen, 
as its values are higher for the odd corrals for the same L, 
and also monotonic with L, unlike even corrals. The spin 
configuration of an open corral can be tuned to a ferrimagnetic 
(FI) {h < he) or planar non-collinear (NC) state {h > he) by 
adjusting N, L or h. (c): closed corrals with even (d): 
closed corrals with odd N. The horizontal axis is N. Due to 
translational invariance the exact value of L does not affect 
the magnetic behavior. For even N the corral changes from a 
FI to a NC configuration with increasing h, while for odd N 
it first goes from a NC to a non-coplanar (NCP) and then to 
a NC configuration. 



corral with STM. 

Up to now, classical spins have been considered in 
Hamiltonian ([1]). However, as mentioned earlier, typical 
values for the corral spins are estimated to be 5^ = 3/2 or 
2 and magnetization curves for Si < 2 are shown in Fig. 
[3l along with the classical result. The quantum mechan- 
ical magnetization approaches the classical curve with 
increasing s^, quite close already for Si = 2. Therefore 
the qualitative features of the classical calculation are ex- 
pected to survive for the actual Si values, in accordance 
with a similar calculation for open chains (37| . 

Closed even corrals. In a closed corral, in contrast 
to open corrals, due to translational invariance of the 
isolated corral the exact linking point L of the extra 
adatom does not affect the magnetic behavior. For an 
isolated even corral there is no frustration originating 
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FIG. 3: Reduced magnetization M/Msat as a function of 
h/hsat for an open corral with = 9 and an extra adatom 
linked at site L = 8. The classical result corresponds to a 
ferrimagnetic (FI) state for lower fields, and a planar non- 
collinear (NC) state for higher fields. The quantum mechan- 
ical magnetization tends to the classical curve as s increases 
from 1/2 towards the typical values 3/2 and 2. The classi- 
cal susceptibility x (upper left inset) has a discontinuity at 
he = 0.120hsat- The polar angles in units of tt (lower right in- 
set) are on the average turning towards the field. For smaller 
fields above the discontinuity spins that originally point along 
the field direction turn away from the field to increase the ex- 
change energy. Nearest neighbor azimuthal angles differ by 
TT, with the two groups of spins having the same azimuthal 
angle indicated by solid (spins at odd sites and extra spin) 
and dashed (spins at even sites) lines respectively. 



in the closed boundary conditions, and nearest neighbor 
spins are antiparallel with no net magnetization in the 
absence of a field (Fig. [B^c)). The extra adatom provides 
an uncompensated spin and M^^o = and it takes a 
finite magnetic field to destroy the FI configuration and 
generate a NC one, similarly to the previous cases. In 
the NC state spins symmetrically placed with respect to 
the linking adatom sl point in the same direction. The 
dependence of he on N is shown in Fig. EJc). As was 
the case before, Mh=o decreases with respect to Msat as 
1/(A^ + 1), thus he decreases with N. 

Closed odd corrals. For odd TV the periodic boundary 
condition introduces a frustrated configuration even for 
an isolated closed corral with no extra adatoms. It is 
not possible for nearest-neighbor spins to be antiparal- 
lel even when = 0, however the net magnetization is 
still zero. The extra adatom provides an uncompensated 
spin and Mh=o = Si^ but the configuration now is NC, 
unlike the FI configurations found before (Fig. Hl^d)). 
This configuration is susceptible to an infinitesimal mag- 
netic field, again unlike all previous cases with Mh=o 7^ 
(Fig. Sl^b)). In it the extra and the linking adatom are al- 
ways antiparallel, with the extra adatom always pointing 
along the direction of the field. Spins which are symmet- 
rically placed with respect to the linking adatom sl have 
the same polar angle (this is true for all magnetic fields 
up to saturation), but they point in opposite directions 




FIG. 4: A closed corral with A" = 9 and an extra adatom. (a) 
Non-coplanar (NCP) configuration of spins for h = 0.171hsat- 
(b) Reduced magnetization M/Msat as a function of h/hsat- 
The classical susceptibility (inset) has two discontinuities at 
h/hsat = 0.160 and 0.861, which divide the planar non- 
collinear (NC) configurations from the NCP. (c) Nearest 
neighbor correlation functions Si -Si+i as a function of h/hsat^ 
where 1' indicates correlation si -sg. They are symmetric with 
respect to the linking adatom sl of the extra spin se (here 
L — 9 and E — 10). For low fields there are correlations that 
get stronger with h. 



in the azimuthal plane (see also Supp. Mat.). In addi- 
tion, nearest neighbor correlations between specific pairs 
of spins increase with /i, and eventually these spins be- 
come antiparallel (Fig. HJ^c)). These pairs are arranged in 
a "dimer" type of configuration, with every second bond 
increasing or decreasing in correlation strength. From 
the ones which are getting stronger with /i, they increase 
in correlation strength as the linking point L of the ex- 
tra adatom is approached. This is counter-intuitive and 
points to the importance of the single extra spin for the 
change of the magnetic properties of the whole corral. 

With increasing field, a susceptibility discontinuity 
leads for the first time to a non-coplanar (NCP) low- 
est spin configuration (Fig. [H^b)), shown in Fig. HJ^a) for 
h = O.lTlhsat- The nearest neighbor correlations now 



decrease only with the magnetic field (Fig. HJ^c)), even 
though there are initially polar angles that do not de- 
crease with /i, due to the competition of exchange and 
magnetic energy (see also Supp. Mat.). For higher fields 
close to saturation a second susceptibility discontinuity 
leads to a NC phase. Now spins symmetrically placed 
with respect to the linking adatom sl point in the same 
direction. The transition fields and the range of exis- 
tence of the different phases as function of N are shown 
in Fig. Efd). The odd closed corral with the extra 
adatom combines the competition between exchange and 
magnetic energy with the frustration introduced by the 
closed boundary conditions, supporting a NCP configu- 
ration not found in any of the other cases. In addition, 
it supports two susceptibility discontinuities in its mag- 
netization curve. Therefore its magnetic behavior is in 
direct contrast with its even counterpart and the open 
corrals as well. 

To conclude, the magnetic properties of corrals (open 
or closed) with an extra adatom attached are differ- 
ent from their counterparts without the extra adatom. 
The extra spin changes the geometrical topology of the 
nanosystem, which has an impact on its magnetic behav- 
ior as a whole. It can tune the lowest energy configura- 
tion of the spins according to the parity of the corral, 
the parity of the linking point of the extra spin and the 
presence or not of periodic boundary conditions. This 
leads to a plethora of configurations that can also be 
three-dimensional, and can be accessed and of use exper- 
imentally. 
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I. SUPPLEMENTAL MATERIAL 

A. Extra Spin Linked with Varying Strength 

The robustness of the infiuence of the extra adatom on 
the isolated corral can be shown by varying the strength 
of the coupling sl-se, which is equal to J in Hamiltonian 
(1) of the main text. If this coupling is scaled with a, 
then the magnetization of an TV = 10 open corral with an 
extra adatom attached at site L = 9 is shown in Fig. [5l 
Even for small values of a the ferrimagnetic state, absent 
when a = 0, appears for small magnetic fields. 

B. Comparison of Classical and Quantum Results 
for Moderate Si 

In Ref. [37] of the main text it has been shown that 
for open isolated corrals without an attached adatom the 
classical result describes accurately the physics up to rel- 
atively small quantum numbers, 5^ > 3/2. As shown in 
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FIG. 5: Reduced magnetization M/Msat as a function of 
h/hsat for an open corral with = 10 and an extra spin 
linked at site L = 9, where the strength J of the sl • se bond 
in Hamiltonian (1) of the main text is scaled with a. The 
ferrimagnetic part of the magnetization curve (low fields) is 
robust, even for small values of a. 
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FIG. 6: Nearest neighbor correlation functions Si ■ s^+i for a 
N = 9 closed corral with an extra adatom attached for Si = 2 
as a function of h/hsat, where l' indicates correlation si • sg. 
The results are similar to the classical results (compare Fig. 
4(c) of the main text). 



Fig. 3 of the main text this is similar for the case of the 
corrals with an extra adatom, as the classical magnetiza- 
tion is closely approached already for such Si values. In 
Fig. 4(c) of the main text the classical nearest neighbor 
correlation functions are plotted as a function of h/hsat- 
The corresponding plot for 5^ = 2 is shown in Fig. [6l 
where by comparison it is seen that the classical result 
again already provides a very good description. 



C. Polar Angles for the = 9 Closed Corral 

The polar angles for the closed corral with N = 9 and 
an extra adatom are shown in Fig. [71 Polar angles sym- 
metrically placed with respect to the linking spin sl are 
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FIG. 7: Polar angles in units of tt as a function of h/hsat 
for the N — 9 closed corral with a single spin attached. The 
two susceptibility discontinuities at h/hsat — 0.160 and 0.861 
(shown with red arrows) separate the three different lowest 
energy configurations, which change from NC to NPC and 
then to NC with increasing field. 



equal. The three different regimes for the lowest energy 
configuration as a function of h/hsat are clearly seen, 
separated by two susceptibility discontinuities. The de- 
pendence of the 9i on h is in general non-monotonic. 



D. Magnetic Anisotropy 



The influence of easy-axis magnetic anisotropy can be 
seen in Fig. [H where a term — ^[X]^i(5f)^ + (-^1;)^] is 
added in Hamiltonian (1) of the main text, with the easy 
axis parallel to the magnetic field axis. For a A/" = 9 
closed corral with an extra adatom and weak anisotropy 
K = 0.01 J the magnetic response does not change sig- 
nificantly. 



J. P. Bucher, D. C. Douglass, and L. A. Bloomfield, Phys. 
Rev. Lett, 66, 3052 (1991). 

J. T. Lau, A. Fohlisch, R. Nietubyc, M. Reif, and 

W. Wurth, Phys. Rev. Lett, 89, 057201 (2002). 

H. Manoharan, C. Lutz, and D. Eigler, Nature 403, 512 

(2000). 

P. Gambardella, A. Dallmeyer, K. Maiti, M. C. Malagoh, 
W. Eberhardt, K. Kern, and C. Carbone, Nature 416, 
301 (2002). 

N. Knorr, H. Prune, M. Epple, A. Hirstein, M. A. Schnei- 
der, and K. Kern, Phys. Rev. B 65, 115420 (2002). 
F. Silly, M. Pivetta, M. Ternes, F. Patthey, J. P. Pelz, 
and W.-D. Schneider, Phys. Rev. Lett 92, 016101 (2004). 

C. F. Hirjibehedin, C. P. Lutz, and A. J. Heinrich, Sci- 
ence 312, 1021 (2006). 

S. Loth, S. Baumann, C. P. Lutz, D. M. Eigler, and A. J. 
Heinrich, Science 335, 196 (2012). 

D. M. Eigler and E. K. Schweizer, Nature 344, 524 
(1990). 

F. J. RueB, W. Pok, T. C. G. Reusch, M. J. Butcher, 
K. E. J. Goh, L. Oberbeck, G. Scappucci, A. R. Hamil- 
ton, and M. Y. Simmons, Smah 3, 563 (2007). 
C. R. Moon, L. S. Mattos, B. K. Foster, G. Zeltzer, and 
H. C. Manoharan, Nat. Nano. 4, 167 (2009). 
F. Meier, L. Zhou, J. Wiebe, and R. Wiesendanger, Sci- 
ence 320, 82 (2008). 

A. J. Heinrich, J. A. Gupta, C. P. Lutz, and D. M. Eigler, 
Science 306, 466 (2004). 

S. Lounis, P. Mavropoulos, P. H. Dederichs, and 

S. Bliigel, Phys. Rev. B 72, 224437 (2005). 

M. S. Ribeiro, G. B. Correa, A. Bergman, L. Nordstrom, 

O. Eriksson, and A. B. Klautau, Phys. Rev. B 83, 014406 

(2011). 

V. S. Stepanyuk, L. Niebergall, W. Hergert, and 

P. Bruno, Phys. Rev. Lett. 94, 187201 (2005). 

A. Bergman, L. Nordstrom, A. B. Klautau, S. Frota- 



XO.16 




FIG. 8: Magnetic susceptibility x for a = 9 closed corral 
with an extra adatom, without and with magnetic anisotropy 
K, with the easy axis parallel to the direction of the mag- 
netic field. A small magnetic anisotropy K = 0.0 IJ does not 
change significantly the magnetic response, which has two 
susceptibility discontinuities at h/hsat — 0.160 and 0.861 for 
K = and at h/hsat = 0.184 and 0.861 for K = O.IJ. 



Pessoa, and O. Eriksson, Phys. Rev. B 73, 174434 (2006). 
[18] J. C. Slonczewski, J. Magn. Magn. Mater. 159, LI 
(1996). 

[19] S. L Kiselev, J. C. Sankey, L N. Krivorotov, N. C. Em- 
ley, R. J. Schoelkopf, R. A. Buhrman, and D. C. Ralph, 
Nature 425, 380 (2003). 

[20] S. Krause, L. Berbil-Bautista, G. Herzog, M. Bode, and 
R. Wiesendanger, Science 317, 1537 (2007). 

[21] M. Menzel, Y. Mokrousov, R. Wieser, J. E. Bickel, 
E. Vedmedenko, S. Bliigel, S. Heinze, K. von Bergmann, 
A. Kubetzka, and R. Wiesendanger, Phys. Rev. Lett. 



7 



108, 197204 (2012). 
[22] M. N. Leuenberger and D. Loss, Nature 410, 789 (2001). 
[23] F. Troiani, A. Ghirri, M. Affronte, S. Carretta, P. San- 

tini, G. Amoretti, S. Piligkos, G. Timco, and R. E. P. 

Winpenny, Phys. Rev. Lett. 94, 207208 (2005). 
[24] A. A. Khajetoorians, J. Wiebe, B. Chilian, and 

R. Wiesendanger, Science 332, 1062 (2011). 
[25] L. Zhou, J. Wiebe, S. Lounis, E. Vedmedenko, F. Meier, 

S. Bliigel, P. H. Dederichs, and R. Wiesendanger, Nature 

Phys. 6, 187 (2010). 
[26] A. A. Khajetoorians, J. Wiebe, B. Chihan, S. Lounis, 

S. Bliigel, and R. Wiesendanger, Nature Phys. 8, 497 

(2012). 

[27] S. Lounis, P. H. Dederichs, and S. Bliigel, Phys. Rev. 

Lett. 101, 107204 (2008). 
[28] S. Holzberger, T. Schuh, S. Bliigel, S. Lounis, and 



W. Wulfhekel, Phys. Rev. Lett. 110, 157206 (2013). 
[29] M. A. Ruderman and C. Kittel, Phys. Rev. 96, 99 (1954). 
[30] T. Kasua, Prog. Theor. Phys. 16, 45 (1956). 
[31] K. Yosida, Phys. Rev. 106, 893 (1957). 
[32] S. Lounis, P. Mavropoulos, P. H. Dederichs, and 

S. Bliigel, Phys. Rev. B 73, 195421 (2006). 
[33] D. Coffey and S. A. Trugman, Phys. Rev. Lett. 69, 176 

(1992). 

[34] N. P. Konstantinidis, Phys. Rev. B 76, 104434 (2007). 
[35] N. P. Konstantinidis, Phys. Rev. B 72, 064453 (2005). 
[36] N. P. Konstantinidis, Phys. Rev. B 80, 134427 (2009). 
[37] A. Machens, N. P. Konstantinidis, O. Waldmann, 

L Schneider, and S. Eggert, Phys. Rev. B 87, 144409 

(2013). 



